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ABSTRACT 


Program  REFINE  was  initially  developed  as  part  of  an  adaptive  finite 
element  modelling  package  for  VAST  where  from  estimates  of  discretization 
errors  based  on  element  residuals,  the  finite  element  model  would  be 
selectively  refined  (by  element  subdivision)  until  acceptable  levels  of 
accuracy  (according  to  computed  error  estimates)  were  attained  [1]. 
However,  increasingly  the  REFINE  program  finds  application  as  a  general 
modelling  tool  which  complements  the  data  generation  capabilities  of 
PVAST  and  HVAST.  Under  these  circumstances,  the  PVAST  and  HVAST  programs 
need  only  concentrate  on  achieving  accurate  geometric  models,  since 
subsequent  refinement  of  the  model  to  attain  acceptable  finite  element 
accuracy  can  be  achieved  with  REFINE.  This  is  particularly  useful  when 
the  initial  model  is  defined  in  a  pointwise  fashion  with  a  graphics 
table.  In  this  context  as  well,  selective  model  refinement  would 
generally  be  employed  but  the  experience  of  the  modeller  would  dictate 
where  the  model  refinements  would  be. 


RESUME 


Le  logiciel  REFINE  a  Ste  conqu  dans  le  cadre  d’un  progiciel  adaptatif  de 
modeiisation  aux  elements  finis  pour  VAST.  Le  programme  raffine 
select ivement  (par  subdivision  des  elements)  le  module  aux  elements  finis 
a  partir  d' estimations  des  erreurs  de  discretization  bashes  sur  les 
r£sidus  des  elements,  jusqu'a  ce  que  le  niveau  de  precision  indique  par 
les  estimations  d'erreur  calcuiees  soit  atteint  [1],  Toutefois,  les 
applications  du  programme  REFINE  se  multiplient  en  tant  qu'outil  general 
de  donn£es  de  PVAST  et  de  HVAST.  Dans  ces  circonstances,  les  programmes 
PVAST  et  HVAST  peuvent  n'etre  utilises  que  pour  produire  des  modeles 
geomdtriques  precis,  REFINE  permettant  de  raffiner  ulterieurement  le 
modele  pour  obtenir  des  elements  finis  de  precision  acceptable.  Cela  est 
particulidrement  utile  lorsque  le  modele  initial  est  d^fini  point  par 
point  a  l’aide  d’une  tablette  graphique.  Dans  ce  contexte  6galement,  le 
moddlisateur  aurait  recours  au  raffinement  seiectif  du  modele,  mais  son 
experience  lui  dicterait  ou  apporter  les  raf f inements . 
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1.1 


CHAPTER  1 
INTRODUCTION 


1 . 1  Introduction 


Program  REFINE  was  initially  developed  as  part  of  an  adaptive  finite 
element  modelling  package  for  VAST  where  from  estimates  of  discretization 
errors  based  on  element  residuals,  the  finite  element  model  would  be 
selectively  refined  (by  element  subdivision)  until  acceptable  levels  of 
accuracy  (according  to  computed  error  estimates)  was  attained  [1]. 
However,  increasingly  the  REFINE  program  finds  application  as  a  general 
modelling  tool  which  complements  the  data  generation  capabilities  of 
PVAST  and  HVAST.  Under  these  circumstances,  the  PVAST  and  HVAST  programs 
need  only  concentrate  on  achieving  accurate  geometric  models,  since  sub¬ 
sequent  refinement  of  the  model  to  attain  acceptable  finite  element 
accuracy  can  be  achieved  with  REFINE.  This  is  particularly  useful  when 
the  initial  model  is  defined  in  a  pointwise  fashion  with  a  graphics 
tablet.  In  this  context  as  well,  selective  model  refinement  would  gener¬ 
ally  be  employed  but  the  experience  of  the  modeller  would  dictate  where 
the  model  refinements  would  be. 

This  report  describes  recent  work  completed  to  extend  the  capabi¬ 
lities  of  the  program  REFINE.  Chapter  2  describes  the  program  modifica¬ 
tions  which  enable  the  refinement  of  additional  element  types.  Almost 
all  element  types  are  operational  now.  Chapter  3  describes  the  intro¬ 
duction  of  appropriate  cubic  multi-point  constraint  equations  for  refine¬ 
ment  of  stiffened  or  unstiffened  panels.  Chapter  4  summarizes  defi¬ 
ciencies  in  refinement  of  substructured  models  and  discusses  how  one  of 
them  was  addressed  by  automatically  defining  newly  generated  boundary 
nodes  (on  superelement  boundaries)  as  either  slave  nodes  or  master  nodes 
in  a  manner  that  preserves  maximum  structural  flexibility.  Chapter  5 
describes  a  restart  capability  for  REFINE  suitable  for  either  substruc— 
tured  or  unsubstructured  finite  element  models.  Chapter  6  outlines  other 
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1 .2 


improvements  in  the  superelement  refinement  algorithm.  Chapter  7  des¬ 
cribes  several  miscellaneous  changes  to  the  REFINE  program  not  called  for 
in  the  contract  statement  of  work  but  which  were  determined  to  be  desir¬ 
able  in  recent  applications  of  the  program. 
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2.1 


CHAPTER  2 

REFINEMENT  OF  ADDITIONAL  ELEMENT  TYPES 


2.1  Introduction 


The  element  refinement  algorithm  of  the  program  REFINE  has  been 
extended  to  several  new  elements  from  the  VAST  element  library  including 
the  curved  beam  element  (IEC=7),  the  axisymmetric  shell  (IEC=19)  and  the 
second  degenerate  form  of  the  brick,  namely  the  square  based  pyramid 
(IEC=2).  In  general,  element  refinement  involves  the  generation  of  new 
nodes,  the  removal  of  the  connectivity  of  the  parent  element  to  be 
refined,  the  generation  of  connectivities  for  new  elements  produced  by 
the  refinement  process  and  also  the  generation  and  manipulation  of  con¬ 
straint  equations  for  irregular  nodes  on  the  interface  between  refined 
elements  and  unrefined  elements.  The  program  modifications  for  each  of 
the  above  mentioned  element  types  will  be  discussed  separately  in  the 
following  sections. 

2.2  Square  Based  Pyramid  (IEC=2) 

The  square  based  pyramid  is  a  special  case  of  the  solid  element  and 
is  defined  by  appropriate  repetition  of  node  numbers  in  the  element 
connectivity  list  (see  Figure  2.1). 

The  generation  of  new  nodes  for  the  square  based  pyramid  followed  an 
approach  similar  to  that  taken  for  the  refinement  of  the  tetrahedron 
elements  and  posed  no  great  difficulties.  The  generation  and  manipula¬ 
tion  of  the  constraint  equations  also  did  not  create  any  unusual  prob¬ 
lems.  The  generation  of  the  connectivities  produced  by  the  element 
refinement  however,  is  a  complex  procedure.  The  philosopy  of  the  REFINE 
program  has  been  that  the  refinement  of  an  element  normally  results  in 
the  generation  of  only  smaller  elements  of  the  same  type.  The  one  excep¬ 
tion  to  this  is  the  refinement  of  the  transition  element  which  results  in 
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2.2 

the  creation  of  transition  and  shell  elements  as  will  be  discussed  in  the 
following  section.  The  refinement  of  the  square  based  pyramid  is  also  an 
exception  to  this  rule  since  the  refinement  results  in  the  formation  of 
square  based  pyramids  and  tetrahedron  elements.  Therefore,  a  new  element 
group  is  necessarily  produced.  The  treatment  of  this  new  group  is 
similar  to  the  new  group  formed  by  refinement  of  the  transition  element. 
The  refinement  of  the  square  based  pyramid  poses  an  additional  problem 
since  the  total  number  of  elements  produced  does  not  follow  the  general 
equation: 

NELMF  =  NOR  IDIEC  (2.1) 

where 

NELMF  =  total  number  of  elements  produced  by  refinement 
NOR  =  order  of  refinement 

IDIEC  =  spatial  dimension  of  the  element  type  (1,2  or  3) 

In  this  case,  IDIEC  is  3  since  the  square  based  pyramid  is  a  three- 
dimensional  element.  The  number  of  pyramid  elements  produced  is  given  by 
the  following  equation: 

NELMP  =  |  *  [N0RIDIEC  +  (2.2) 

where 

NELMP  =  total  number  of  pyramid  elements  produced  by  refinement 
NOR  =  order  of  refinement 
IDIEC  =  3 

The  number  of  tetrahedron  elements  produced  is  given  by  the  following 
equation: 


NELMT  =  |  [N0RIDIEC 


-  NOR] 


(2.3) 


P151394.PDF  [Page:  14  of  80] 


2.3 


where 

NELMT  =  total  number  of  tetrahedron  elements  produced  by  refinement 
NOR  =  order  of  refinement 
IDIEC  =  3 

The  refinement  of  the  square  based  pyramid  does  produce  an  equivalent 
number  of  elements  to  Equation  2.1  if  half  the  number  of  tetrahedrons 
given  by  Equation  2.3  are  added  to  the  total  number  of  pyramids  given  by 
Equation  2.2.  This  manipulation  is  somewhat  justified  with  the  observa¬ 
tion  that  each  pyramid  may  be  broken  into  2  component  tetrahedrons. 
Table  1  summarizes  element  information  for  several  different  orders  of 
refinement.  It  can  be  readily  seen  that  combining  Equation  2.2  with  half 
of  Equation  2.3  always  yields  Equation  2.1. 

Now  that  the  number  and  type  of  elements  produced  have  been  estab¬ 
lished,  the  next  consideration  is  the  orientation  of  the  elements.  The 
orientation  of  elements  is  important  when  generating  the  load  data  for 
the  refined  model.  The  new  pyramid  elements  are  therefore  generated  with 
their  face  orientation  the  same  as  the  parent  element.  The  new  tetra¬ 
hedron  elements  produced  on  each  of  the  faces  have  consistently  one 
orientation.  The  logic  required  in  the  program  to  generate  the  element 
connectivities  is  quite  complex  and  being  sufficiently  unique,  it  forms  a 
separate  section  in  subroutine  REFINE.  Load  refinement  for  square  based 
pyramids  has  not  been  performed  under  this  contract. 

2.3  Transition  Element  (IEC=6) 


Refinement  of  the  transition  element  presents  some  unique  difficul¬ 
ties  and  some  further  development  is  still  required.  Presently,  a 
transition  element  refinement  of  degree  n  generates  n  transition  elements 
and  (n-l)xn  shell  elements.  The  total  number  of  elements  produced  is 
governed  by  Equation  2.1.  As  illustrated  in  Figure  2.2  for  a  refinement 
degree  of  2,  the  face  adjoining  a  solid  element  is  divided  into  n  ele- 
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ments  with  no  element  divisions  through  the  element  thickness.  Element 
refinements  of  the  adjacent  solid  element  with  refinement  degree  n,  on 
the  otherhand,  results  in  n  element  divisions  through  the  thickness  and 
thus  a  refinement  inconsistency  on  the  common  interface.  In  fact, 
depending  whether  the  transition  or  solid  element  is  refined  first, 
difficulties  are  sometimes  experienced  in  generating  correct  multi-point 
constraint  equations  on  the  transition-to-solid  element  interface. 

Although  a  coding  error  was  previously  suspected,  as  the  cause  of 
the  difficulty  in  generating  correct  MPC  equations,  this  has  now  been 
ruled  out.  After  considerable  study,  it  became  apparent  that  for  reli¬ 
able  MPC  equations  the  MPC  on  the  transition-to-solid  element  interface 
modification  of  the  MPC  equation  generation  algorithm  for  the  solid 
element  was  unavoidable.  Different  MPC  equations  would  be  generated  on 
the  faces  of  a  solid  element  depending  on  whether  a  transition  element  or 
another  solid  element  attached  to  it.  In  order  to  achieve  this  requires 
considerable  changes  to  the  databases  on  which  REFINE  operates.  Infor¬ 
mation  related  to  the  neighbours  of  each  element  in  the  model  would  have 
to  be  generated  and  manipulated  during  refinement  of  a  model.  It  is  not 
difficult  to  appreciate  the  difficulty  of  the  task  especially  when  one 
contemplates  the  refinement  of  a  model  which  has  been  refined  previously 
and  includes  MPC  equations  on  the  interface  between  refined  and  unrefined 
sections.  The  philosopy  of  the  REFINE  program  has  been  that  each  element 
selected  by  the  user  can  be  refined  and  appropriate  MPC  equations 
generated  without  regard  to  its  neighbours. 

An  alternative  refinement  strategy  involving  the  division  of  a 
transition  element  into  smaller  transition  elements  and  wedge  elements  as 
illustrated  in  Figure  2.3  has  therefore  been  considered.  Although  this 
has  the  advantages  that  any  given  interface  is  always  divided  consistent¬ 
ly  and  that  the  solid  element  zone  moves  slightly  into  the  shell  struc¬ 
ture,  some  deterioration  of  element  proportioning  develops  if  extensive 
refinements  in  the  original  transition  element  are  required.  However,  if 
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the  maximum  refinement  order  is  no  more  than  five  or  six  and  the  maximum 
aspect  ratio  (element  length  to  thickness)  is  no  greater  than  ten,  e 
aspect  ratio  will  be  acceptable. 


The  work  performed  under  this  contract  consisted  of  a  review  of  the 
existing  coding  to  determine  the  most  appropriate  treatment  of  the 
transition  element.  It  can  be  concluded  that  refinement  of  transition 
elements  cannot  be  provided  within  the  framework  of  REFINE  without 
compromising  the  assumptions  on  which  the  program  presently  operates. 
The  first  refinement  strategy  involving  changes  to  the  program  data 
structures  represents  a  major  undertaking  and  should  not  be  adopted 
unless  adequate  manpower  can  be  budgeted. 

2.4  Curved  Beam  Element  (IEG*7)  and  Axisymmetric  Shell  (IECf19I 


The  curved  beam  element  (superparametric  beam  element)  (IEC=7)  is 
suitable  as  a  stiffener  for  the  thick/thin  shell  element  <IEC=1)  and  if 
used  as  such,  the  displacement  node  numbers  of  the  shell  to  which  the 
beam  is  connected  are  to  be  identified.  These  nodes  must  be  identified 
in  the  REFINE  program  since  any  constraint  equations  generated  must 
relate  to  the  displacement  nodes  and  not  to  the  geometric  nodes, 
same  procedure  for  identifying  the  displacement  nodes  is  also  necessary 
for  the  axisymmetric  shell  element  (IEC=19>  since  the  constraint  equa¬ 
tions  relate  to  those  nodes.  The  refinement  of  the  curved  beam  elements 
also  required  the  development  of  a  new  subroutine  CBEAM  to  calculate  the 
curved  beam  properties  at  new  nodes  using  the  shape  function  of  each  new 

nodal  point  in  question. 
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TABLE  2.1 

Number  of  Elements  Produced  by 
Refinement  of  a  Square  Based  Pyramid 


NOR 

Order  of 
Refinement 

NELMF 

General  IDIEC=3 
(Equation  2.1) 

NELMP 

Pyramid  Elements 
(Equation  2.2) 

NELMT 

Tetrahedron  Elements 
(Equation  2.3) 

1 

1 

1 

0 

2 

8 

6 

4 

3 

27 

19 

16 

4 

64 

44 

40 
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CUBE  MODELLED  BY  SQUARE  BASED  PYRAMID  ELEMENTS 
O 


27  27  I 


1 

10.0000 

0.0000 

10.0000 

2 

10.0000 

0.0000 

5. OOOO 

o 

o 

10.0000 

0.0000 

0.0000 

4 

10.0000 

5. 0000 

10. oooo 

5 

10.0000 

5. 0000 

5. 0000 

6 

10.0000 
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FIGURE  2.1:  Typical  GOM  File  for  VAST  With  Square-Based  Pyramid 
Elements . 
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CHAPTER  3 

CUBIC  MULTIPOINT  CONSTRAINT  EQUATIONS  FOR  REFINEMENT 
OF  STIFFENED  OR  UNSTIFFENED  PANELS  BY  PROGRAM  REFINE 


3.1  Introduction 


If  a  finite  element  model  is  refined  nonuniformly,  more  nodes  may  be 
generated  on  one  side  of  some  element  interfaces  within  the  original 
model  than  on  the  other  side.  Nodes  not  joined  to  nodes  in  contiguous 
elements  across  such  an  interface  are  referred  to  as  irregular  nodes. 
The  presence  of  irregular  nodes  on  an  element  interface  implies  lack  of 
displacement  compatibility  unless  displacements  at  all  irregular  nodes 
are  related  to  displacements  at  regular  nodes  on  such  an  interface. 

Constraint  equations  are  typically  of  the  form: 


?t 

udm  ~  Cmn  uin 


where  U(jm  is  the  n-th  dependent  freedom  (either  translational  displace¬ 
ment  or  rotation  component),  Uin  are  the  independent  freedoms  on  which 
they  depend,  and  c^  are  coefficients  relating  them.  For  almost  all 
elements  in  the  VAST  library,  degree  of  freedom  i  at  the  dependent  node  m 
is  related  to  only  degree  of  freedom  i  at  two  or  more  independent  nodes 
n.  What  this  means  for  a  refinement  program  such  as  REFINE  is  that 
constraint  equation  coefficients  need  be  stored  for  only  the  first 
degree-of-freedom  at  each  dependent  node.  At  the  end  of  the  program 
REFINE  when  multi-point  constraint  data  is  being  written  to  the  USE  file 
for  the  refined  model,  constraint  equations  for  all  degrees  of  freedom  at 
dependent  nodes  can  be  created  by  simple  repetition.  The  plate  bending 
element  and  the  general  beam  stiffener  element  however  are  exceptions. 
Out-of-plane  displacements  for  points  along  an  edge  are  cubic  functions 
in  a  linear  edge  variable  (defined  by  values  of  out-of-plane  displacement 
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and  rotations  at  end  points ) .  The  present  note  outlines  the  required 
modifications  to  generalize  constraint  generation  algorithms  within 
program  REFINE. 

A  limited  capability  for  the  generation  of  cubic  multi-point 
constraint  equations  for  the  triangular  plate  bending  element  and  the 
general  beam  stiffener  elements  has  now  been  implemented  and  tested  in 
program  REFINE.  It  permits  nonuniform  refinements  of  stiffened  or 
unstiffened  plate  panels.  The  panel  must  be  flat  and  its  normal  aligned 
with  one  of  the  coordinate  axes.  The  user  is  prompted  for  the  orienta¬ 
tion  of  the  normal.  Although  it  is  relatively  straightforward  to  compute 
the  direction  of  the  normal  for  a  stiffened  panel,  the  program  nonethe¬ 
less  prompts  the  user  to  manually  define  its  direction  cosines.  However, 
it  is  felt  that  in  facetted  plate  structures,  a  single  normal  may  have  to 
be  utilized  for  correct  generation  and  manipulations  of  the  multi-point 
constraint  equations. 

What  distinguishes  constraint  equation  generation  for  the  plate 
bending  and  general  beam  stiffener  element  types  from  others  in  the  VAST 
library  is  that  typically  a  different  constraint  equation  is  required  on 
each  freedom  at  irregular  nodes  of  a  nonuniformly  refined  model. 
Development  of  this  capability  thus  involved  two  principal  steps.  The 
first  step  consisted  of  the  generalization  of  procedures  with  REFINE  so 
that  more  than  one  constraint  equation  could  be  generated  and  manipulated 
for  each  irregular  (constrained)  node.  The  second  step,  of  course,  was 
to  implement  the  appropriate  cubic  constraint  equations  required  for  out- 
of -plane  displacements  in  the  plate  bending  elements.  Similar  approach 
is  necessary  for  the  general  beam  stiffener  element. 

In  the  following  sections  of  this  note  these  steps  are  elaborated  on 
and  some  program  usage  restriction  and  guidelines  are  also  discussed. 
The  future  development  of  this  work  is  also  discussed  briefly. 
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3.2  Generalization  of  Constraint  Equation  Manipulation 

The  basic  logic  implemented  within  program  REFINE  for  generation, 
manipulation,  and  optional  elimination  of  constraint  equations  (as  well 
as  for  computation  of  new  nodal  coordinates,  and  for  generation  of 
connectivities  for  new  elements)  is  adequate  for  processing  of  cubic 
constraint  equations.  The  same  nodes  are  related  in  each  constraint 
equation  associated  with  a  given  irregular  node;  different  degrees  of 
freedom  are  generally  involved  in  each  constraint  equation  and  frequently 
more  than  one  degree  of  freedom  is  used  at  every  dependent  node.  Depen¬ 
dent  (irregular)  nodes  for  which  displacements  are  defined  via  constraint 
equations  are  stored  in  array  NK  and  associated  independent  nodes  in  the 
matrix  NJ.  Enough  information  can  be  extracted  from  the  structure  of  the 
first  constraint  equation  at  every  dependent  node  for  the  purpose  of 
filling  arrays  NK  and  NJ. 

It  is  primarily  in  the  generation  and  handling  of  the  coefficients 
for  the  constraint  equation  that  implementation  of  cubic  constraint  equa¬ 
tions  for  the  plate  bending  elements  and  general  beam  stiffener  element 
differs  from  the  implementation  of  appropriate  constraint  equations  for 
other  elements  in  the  VAST  library.  In  this  section,  the  question  of 
modifications  to  the  constraint  equation  handling  is  discussed. 

The  simplest  method  of  extending  the  program  logic  to  handle 
multiple  constraint  equations  where  conceivably  a  different  constraint 
equation  is  associated  with  each  degree  of  freedom  at  a  dependent  node, 
is  to  leave  the  basic  program  logic  related  to  determining  the  status  of 
nodes  (regular  or  irregular),  to  keeping  track  of  all  independent  nodes 
associated  with  each  dependent  one  etc.  intact  and  setup  an  auxilliary 
mechanism  (data  management  system)  for  storage  and  retrieval  of 
constraint  equation  coefficient  data.  Whenever  constraint  equation 
manipulations  are  required,  instead  of  one  set  of  equations  being 
operated  on  several  sets  will  be.  Generally  the  number  of  constraint 
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equations  to  be  operated  on  will  be  equal  to  the  number  of  nodal  degrees 
of  freedom  at  each  irregular  node.  Of  course,  repeating  constraint 
equation  computations  when  all  constraint  equations  associated  with  an 
irregular  node  are  equivalent  is  unnecessary  and  future  refinements  of 
this  algorithm  could  conceivably  involve  investigating  the  possibility  of 
performing  one  set  of  constraint  equation  calculations  when  they  are 
equivalent  for  each  degree  of  freedom  at  irregular  nodes.  For  such  nodes 
of  course,  it  will  be  necessary  to  expand  out  the  constraint  equations 
and  generate  the  NDF  constraint  equations  from  the  single  set  of 
constraint  equation  coefficients  generated  by  the  program. 

Due  to  efforts  to  improve  speed  of  program  REFINE  during  an  earlier 
contract  [  1  ] ,  dependent  node  numbers  and  associated  constraint  equation 
coefficients  were  stored  in  program  arrays  to  reduce  the  amount  of 
program  I/O  and  thereby  increase  computational  efficiency  of  critical 
sections  of  the  refinement  algorithm.  To  store  coefficients  for  all 
constraint  equations  at  each  irregular  (dependent)  node  will  raise 
program  storage  requirements  quite  dramatically  particularly  since  the 
maximum  number  of  structural  nodes  has  recently  been  increased  to  2500. 
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CHAPTER  4 

AUTOMATIC  TREATMENT  OF  BOUNDARY  NODES  IN  SUPERELEMENT  REFINEMENT 

4.1  Introduction 


Substructuring  plays  a  major  role  in  finite  element  modelling  of 
ship  structures,  and  any  finite  element  model  refinement  capability  would 
be  incomplete  if  it  could  not  handle  substructured  finite  element 
models.  Earlier  versions  of  the  program  REFINE  were  not  fully 
operational  on  superelements  [1].  A  number  of  restrictions  were  adopted 
to  facilitate  producing  a  preliminary  capability  to  refine  superelements 
which  included: 

a.  Only  level  one  superelements  can  be  refined. 

b.  A  substructure  may  only  define  one  superelement . 

c.  Substructures  had  to  be  refined  in  ascending  order. 

d.  All  new  superelement  interface  nodes  had  to  become  either 
master  nodes  or  slave  nodes  as  defined  by  an  interactive 
prompt . 

e.  Limitations  on  refinements  of  previously  refined  structures  are 
permitted,  notably,  each  element  of  the  original  model  could  be 
refined  only  once. 

Considerable  effort  in  this  contract  has  been  concentrated  on  removing 
the  above  mentioned  restrictions.  This  chapter  describes  work  on  Item  d. 

4.2  Importance  of  Automatic  Classification  of  Boundary  Nodes 

Geometrically,  refinement  of  substructured  models  is  no  different 
than  refinement  of  unsubstructured  models.  Basically  element  refinement 
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still  involves  the  generation  of  new  nodes,  removal  of  the  connectivity 
for  the  parent  element,  the  generation  of  connectivities  for  new  elements 
produced  by  element  refinement  as  well  as  generation  and  manipulations  of 
the  constraint  equations  for  irregular  nodes  on  the  interface  between 
refined  and  unrefined  elements. 

Some  auxiliary  issues  related  to  substructures  alone  cannot  be 
avoided  when  they  are  refined.  Firstly,  constraint  equations  are  imple¬ 
mented  by  the  slave  node  concept  whereby  dependent  nodes  for  constraint 
equations  are  master  nodes.  Consequently,  nodes  internal  to  substruc¬ 
tures  may  have  to  become  master  nodes  in  order  that  constraint  equations 
can  be  applied.  Secondly,  assignment  of  constraint  equations  for  bound¬ 
ary  nodes  in  one  substructure  (superelement)  cannot  be  done  independently 
of  other  substructures  of  the  model.  If  all  substructures  sharing  a 
given  interface  are  refined  consistently  all  new  nodes  generated  on  the 
interface  can  become  master  nodes.  Otherwise,  those  new  nodes  which  do 
not  have  counterparts  in  all  connecting  substructures  must  be  slaved 
(have  constraint  equations  retained  for  them).  A  limited  capability  to 
refine  substructure/superelements  was  implemented  in  the  previous 
contract  [1].  The  approach  adopted  was  to  prompt  the  user  whether  all 
boundary  nodes  were  to  be  identified  slave  or  master  nodes.  The  user  was 
thus  responsible  for  ensuring  that  the  refinement  of  all  substructures/ 
superelements  were  performed  appropriately  for  the  selection  made. 
Figure  4.1  illustrates  the  two  modelling  options  for  new  superelement 
interface  nodes. 

The  restrictions  adopted  for  the  initial  substructure/superelement 
refinement  capability  developed  partly  from  time  and  budgetary  con¬ 
straints.  At  the  time,  it  also  seemed  that  the  implementation  of  a 
general  approach  for  automatically  identifying  a  boundary  node  to  be  a 
master  node  or  slave  node  would  require  searching  all  element  connectivi¬ 
ties  of  all  superelements  within  the  structure.  The  large  amounts  of 
file  searching  implied  significant  increases  in  CPU  requirements  and 
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decreased  algorithm  (program)  speed  when  substructuring  was  involved.  In 
fact,  recent  developments  in  subroutine  MPC  lead  to  the  proposal  that  the 
approach  taken  to  determine  when  a  constrained  node  can  be  eliminated 
could  also  be  used  to  determine  when  a  constrained  boundary  node  can  be 
eliminated.  If  the  constraints  on  a  boundary  node  remained,  the  node 
would  be  identified  as  a  slave  node,  and  if  the  constraints  are  removed, 
the  node  would  be  identified  as  a  master  node. 

The  main  procedures  employed  in  subroutine  MPC  to  determine  whether 
constraints  on  nodes  could  be  eliminated  required  a  search  through  all 
the  element  connectivities.  The  extension  of  this  procedure  to  a 
substructured  model  would  thus  require  a  search  through  all  the  element 
connectivities  of  all  the  superelements.  At  the  start  of  this  contract, 
a  review  of  the  data  requirements  for  determining  whether  a  boundary  node 
should  be  specified  as  a  slave  or  master  node  revealed  that  sufficient 
data  for  a  reliable  decision  could,  in  fact,  be  obtained  from  the  super¬ 
element  data  section  of  the  input  file.  This  data  provides  the  master 
node  list  for  each  superelement  and  thus  defines  the  connectivity  of  each 
superelement.  A  search  through  this  data  would  eliminate  the  need  for 
searching  through  all  the  individual  elements  in  the  structure.  The 
implementation  of  this  procedure  within  a  new  subroutine  called  MPCSE  is 
detailed  in  this  chapter.  The  overall  structure  of  the  REFINE  program 
does  not  change  beyond  the  addition  of  a  new  module  as  is  readily  appre¬ 
ciated  by  examining  Figure  4.2. 

The  capability  to  refine  substructure/superelements  has  been 
improved  in  this  contract  to  automatically  identify  the  boundary  nodes  as 
slave  or  master.  A  boundary  node  is  defined  as  a  node  located  on  the 
perimeter  of  a  superelement. 

This  chapter  will  review  the  steps  taken  to  identify  nodes  requiring 
multipoint  constraint  data.  Constraint  equations  for  substructures  are 
implemented  by  using  the  constrained  substructure  option  of  VAST  by 
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defining  independent  nodes  as  master  nodes  and  dependent  nodes  as  slave 
nodes.  The  data  for  constrained  node  is  manipulated  on  scratch  files 
using  the  stiffness  modification  format  and  is  re-formatted  to  the  super¬ 
element  format  only  in  the  final  step.  From  the  point  of  view  of  geome¬ 
try  an  individual  substructure  is  the  same  as  an  unsubstructured  model. 
The  treatment  of  the  multipoint  constraint  data  as  outlined  above  allows 
the  REFINE  program  to  operate  basically  the  same  way  for  both  a  substruc¬ 
tured  and  an  unsubstructured  model.  The  only  difference  is  that  for  a 
substructured  model  the  constraint  equations  are  again  reviewed  after  all 
the  substructures  have  been  refined. 


4.3  Review  of  Important  Terminology 

The  following  definitions  will  be  useful  when  discussing  constraint 
equation  generation  and  manipulation.  The  definitions  assigned  will  be 


used  throughout  this  chapter. 

Border  Node  A  node  which  is 

--  parent  element . 


located  on  the  perimeter  of  a 


Boundary  Node 


A  node  which  is  located  in  the  perimeter  of  a 
superelement . 


Dependent  Node 
Interface  Node 
Internal  Node 


A  node  requiring  a  multipoint  constraint  equation. 

A  border  node  common  to  two  or  more  elements . 

A  node  which  is  not  located  in  the  perimeter  of  a 
superelement . 


4.4  identification  of  Node  Requiring  Multipoint  Constraint  Equations 

Multipoint  constraint  equations  are  required  for  irregular  nodes 
created  by  non-uniform  mesh  refinement.  The  REFINE  program  generates  the 
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constraint  equations  for  all  "border"  nodes  created  as  the  result  of  an 
element  refinement  and  then  may  remove  the  constraint  through  a  series  of 
checks  performed  at  three  stages  before  the  generation  of  the  final  out¬ 
put  file.  The  constrained  node  is  investigated  to  ensure  the  constraint 
is  required  by  subroutine  REFINE,  subroutine  MPC,  and  if  the  model  is 
substructured,  subroutine  MPCSE.  The  REFINE  subroutine  is  called  as  each 
individual  element  is  refined.  Subroutine  MPC  is  called  after  all  the 
elements  in  a  substructure  have  been  refined  and  MPCSE  is  called  after 
all  the  substructures  have  been  refined.  The  checks  performed  by  each  of 
the  subroutines  are  described  in  detail  in  the  following  sections. 

4.5  Constraint  Equation  Operations  in  Subroutine  REFINE 

Subroutine  REFINE  generates  the  new  element  connectivities  for  the 
elements  created  by  the  refinement  of  an  element  and  also  generates  and 
investigates  the  multipoint  constraint  equations.  The  constraint  equa¬ 
tions  are  calculated  by  subroutines  ELEML,  ELEMQ ,  ELEMC  or  ELEMP  depen¬ 
ding  upon  the  element  type.  Subroutine  REFINE  does  not  generate  or 
investigate  those  nodes  that  are  completely  internal  to  the  parent  ele¬ 
ment  and  generates  constraints  only  for  new  nodes  (not  previously  in  the 
nodal  coordinate  array).  However,  if  a  node  is  a  border  node  and  is  not 
new  it  is  checked  to  see  if  it  is  already  specified  as  a  dependent  node. 
If  it  is  a  dependent  node  then  the  equation  is  checked  to  ensure  all  the 
nodes  are  included  in  the  connectivity  of  the  element  being  refined.  If 
the  nodes  are  included  then  the  equation  is  investigated  further  to 
determine  if  the  element  being  refined  is  the  only  element  defining  the 
independent  basis  (for  the  constraint  equation).  This  operation  is  per¬ 
formed  by  calling  subroutine  ELEMR  with  the  ISW  parameter  set  to  three. 
The  constraint  equation  is  checked  against  all  the  original  element  con¬ 
nectivities.  If  the  element  being  refined  is  the  only  element  defining 
the  independent  basis  then  the  constraint  is  removed.  Otherwise  the  node 
is  an  "interface"  node  and  this  is  indicated  by  filling  the  second  posi¬ 
tion  of  the  NIDO  array  with  the  number  of  elements  containing  the  inde— 
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pendent  basis  in  their  connect ivitity.  The  identification  of  the  nodes 
as  an  "interface"  node  is  required  for  subroutine  MPC. 

4.6  Constraint  Equation  Operations  in  Subroutine  MPC 

Subroutine  MPC  determines  which  multipoint  constraint  equations  are 
required  and  eliminates  those  which  are  not.  The  procedure  to  determine 
whether  a  constraint  can  be  removed  is  described  as  follows: 

1.  All  the  elements  connectivities  (if  superelements  are  present,  then 
the  connectivities  for  substructure  being  refined)  are  read.  After 
each  connectivity  is  read,  the  independent  nodes  of  each  multipoint 
constraint  equation  are  compared  against  the  connectivity  nodes  for 
all  elements.  If  all  independent  nodes  are  present  in  any  connecti¬ 
vity  list,  the  equation  must  be  retained.  Otherwise,  further  inves¬ 
tigation  of  the  multipoint  constraint  equation  is  necessary  before 

it  can  be  removed. 

2.  If  all  independent  nodes  of  any  multipoint  constraint  equation  are 
not  contained  in  any  element  connectivity,  the  equation  is  checked 
to  ensure  its  independent  nodes  are  not  in  fact  dependent  (as  a 
result  of  constraint  equations  considered  previously).  If  dependent 
nodes  are  found  and  the  node  is  not  an  "interface"  node  then  the 
constraint  equation  remains.  If  not,  and  the  geometric  data  is  also 
not  substructured,  the  constraint  is  marked  to  be  removed. 

3.  This  step  is  performed  only  if  the  file  is  sub  structured.  The 
independent  nodes  are  checked  to  see  if  they  correspond  to  master 
nodes  and  if  the  master  nodes  are  found  only  in  the  substructure 
being  refined,  then  the  constraint  is  removed.  If  the  nodes  do  not 
correspond  to  master  nodes  and  the  node  is  marked  as  an  interface 
node  then  the  constraint  is  also  removed.  This  step  also  serves  to 
identify  "boundary"  nodes  which  are  required  for  subroutine  MPCSE 
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since  the  constraint  may  still  yet  be  removed  depending  upon  the 
refinement  of  other  substructures. 

After  the  above  steps  are  performed  for  an  unsubstructured  model, 
the  constraints  which  remain  are  condensed  and  skew  effects  are  taken  in¬ 
to  account.  The  resulting  final  data  is  written  to  the  output  SMD  file. 

For  a  substructured  file,  additional  steps  are  performed  and  the 
file  data  is  not  condensed.  The  "boundary"  nodes  identified  in  step  3 
have  their  coordinate  locations  compared  to  the  existing  list  of  boundary 
nodes  to  determine  their  boundary  node  number.  This  procedure  is  per¬ 
formed  by  subroutine  IDBLIS.  The  master  node  array  for  the  substructure 
being  refined  is  then  filled  with  the  negative  "boundary"  node  number. 
Next,  the  independent  nodes  are  compared  to  the  master  node  list  and  if 
they  are  not  present  then  the  list  is  incremented.  In  the  case  of  a 
boundary  node,  the  boundary  node  number  and  the  master  node  number  it  is 
dependent  upon  are  written  to  the  NBOR  file  for  processing  in  subroutine 

MPCSE.  The  data  for  all  the  constrained  nodes  is  written  to  the  NREB 
file. 

^•7  Constraint  Equation  Operations  in  Subroutine  MPCSE 

After  the  refinement  of  all  elements  in  all  substructures  has  been 
completed,  subroutine  MPCSE  is  called.  At  this  point,  all  the  multipoint 
constraint  equations  have  been  computed  and  stored  on  the  NREB  file. 
These  equations  are  adjusted  to  take  into  effect  skew  conditions  but  are 
not  condensed.  This  subroutine  makes  use  of  existing  array  space.  The 
steps  to  determine  whether  the  constraint  is  to  remain  are  as  follows: 

1.  A  loop  over  the  superelements  is  performed  to  read  the  data  on  the 
NBOR  file  and  fill  the  NJ  array  with  the  master  node  numbers  that 
the  boundary  node  is  dependent  upon.  The  NBOR  file  was  generated  in 
subroutine  MPC  and  has  the  format  given  in  Table  4.1.  The  substruc- 


ture  and  master  node  numbers  for  each  substructure  are  written  to  a 
file  and  a  second  file  containing  the  boundary  node  numbers  is 

.created . 


A  loop  over  the  superelements  is  again  performed  to  read  the  master 
node  numbers  in  core  and  compare  them  against  the  array  generated  in 
step  1  to  determine  if  the  superelement  contains  all  the  master 
nodes  that  the  boundary  node  depends  upon.  If  YES,  then  another 
array,  IKEC,  is  filled  with  the  superelement  number  and  the  counter 
for  the  number  of  superelements  containing  the  independent  basis  of 
the  equation  is  incremented. 

The  superelements  are  again  looped  over  and  the  boundary  nodes  are 
now  placed  in  core  and  compared  against  the  array  generated  in  step 
2  The  constraint  remains  on  a  node  if  the  independent  basis  is 
contained  in  the  superelement  as  flagged  in  the  1REC  array  but  the 
superelement  does  not  contain  the  boundary  node  as  determined  by  the 
superelements  boundary  node  list.  The  one  exception  is  if  the 
independent  basis  is  found  only  in  the  superelement  then  it  is 
further  investigated  to  determine  if  the  node  is  an  interior  node. 
If  the  node  is  an  interior  node,  it  is  neither  a  slave  or  master 

node. 

New  master  node  numbers  are  assigned  to  all  the  border  nodes  flagged 
as  having  their  constraint  removed. 


The  superelements  are  again  looped  over  and  the  following  steps  are 
performed: 

(a)  The  boundary  node  numbers  flagged  as  negatives  in  the  master 
node  list  for  each  superelement  are  substituted  with  the  new 
master  node  number  calculated  in  Step  4. 
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(b)  The  multipoint  constraint  data  generated  by  subroutine  MPC  is 
read  and  the  data  is  checked  to  see  which  constraints  may  be 
removed.  If  the  node  was  determined  to  be  a  master  node,  then 
the  constraint  is  removed  at  this  point. 

(c)  The  multipoint  constraint  equations  to  remain  are  condensed. 

(d)  The  data  is  transformed  from  the  ISTIFM  data  format  to  the 
IELEMB  format  and  written  to  the  output  file. 

4.8  Sample  Problems 

In  this  section,  the  results  of  refining  an  unsubstructured  model  is 
first  presented  to  demonstrate  how  the  REFINE  program  determines  which 
nodes  to  have  multi-point  constraint  equations  applied.  A  substructured 
model  is  then  refined  to  illustrate  the  treatment  of  boundary  nodes. 

Unsubstructured  Model  Prob lem 


The  unrefined  model  is  a  square  plate  composed  of  nine  4— noded  quad¬ 
rilateral  shell  elements  (IEC=5)  arranged  in  a  3x3  mesh.  The  central 
element  is  first  refined  to  order  2  and  then  two  of  the  resulting 
elements  are  refined  again  to  order  2.  The  resulting  mesh  is  shown  in 
Figure  4.3  and  the  following  discussion  will  refer  to  it. 

The  refinement  of  the  central  element  results  in  new  nodes  17 
through  21.  The  border  nodes  17,  18,  20  and  21  are  specified  as  con¬ 
strained  in  subroutine  REFINE. 

The  refinement  of  the  first  new  element  (with  corner  nodes  6,  17, 
19,  18)  results  in  new  nodes  22  through  26.  The  border  nodes  22,  23,  25 
and  26  are  specified  as  constrained  in  subroutine  REFINE. 
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The  refinement  of  the  second  new  element  (with  corner  nodes  17,  7, 
20,  19)  results  in  new  nodes  27  through  30.  The  border  nodes  27,  29  and 
30 ?  are  specified  as  constrained  in  Subroutine  REFINE.  Node  25  is  a  bor¬ 
der  node  but  is  not  new  and  is  determined  to  have  been  already  specified 
as  a  constrained  node.  The  node  is  specified  as  an  -interface-  node 
since  its  independent  basis  is  contained  in  more  than  one  element. 

The  element  refinement  specifications  are  complete  and  now  subrou¬ 
tine  MPC  is  used  to  determine  which  constraint  equations  are  required. 
The  constraint  equations  generated  are  uncondensed  which  means  a  con¬ 
strained  node  is  related  to  its  two  adjacent  nodes  along  the  edge  being 
refined.  Table  4.2  identifies  the  independent  nodes  for  each  dependent 
node  for  the  uncondensed  equations  of  this  sample  problem.  Constrained 
nodes  17,  18,  20  and  21  remain  since  the  comparison  against  the  connecti¬ 
vity  list  revealed  that  the  independent  nodal  basis  for  each  is  contained 
in  the  connectivity  of  at  least  one  element  in  the  model.  Constrained 
nodes  22,  23,  26,  27,  29  and  30  also  remain  since  the  independent  nodal 
basis  of  each  contains  dependent  nodes  which  are  not  interface  nodes. 
The  constraint  equations  on  node  25  were  removed  at  this  point  since  the 
node  was  flagged  as  an  interface  node. 

The  constrained  nodes  to  remain  are  condensed  and  the  resulting  data 
is  written  to  the  input  file. 

Substructured  Model  Problem 

The  unrefined  model  is  a  ship  transverse  bulkhead  composed  of  four 
structures/superelements  shown  in  Figure  4.4.  Both  superelements  1  and  2 
will  be  partially  refined  along  their  interface  and  the  resulting  mesh  is 
shown  in  Figure  4.5.  VASTG  graphics  has  no  provision  for  identification 
of  constrained  nodes  (slave  nodes)  for  substructures.  It  is  therefore 
not  possible  to  demonstrate  the  operation  of  REFINE  for  substructures 
using  graphics  only  and  necessary  reference  will  be  made  to  the  super¬ 
element  data  ( PREFX . SED ) . 
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The  superelement  data  for  the  unrefined  model  is  shown  in  Table 
4.2.  Superelement  1  has  64  master  nodes  which  are  plotted  in  Figure  4.6 
using  the  option  to  plot  border  nodes  labelled  with  local  node  numbers 
and  superelement  2  has  80  nodes  which  are  plotted  in  Figure  4.7-  The 
refinement  of  superelements  1  and  2  result  in  the  generation  of  models 
with  nodes  as  shown  in  Figures  4.8  and  4.9.  The  new  nodes  generated 
internally  to  the  refined  mesh  create  no  changes  to  the  number  of  slave 
or  master  nodes.  Table  4.5  presents  the  refined  superelement  data  for 
the  model.  The  ten  nodes  in  each  refined  superelement  which  are  identi¬ 
fied  as  slave  nodes  in  this  data  have  been  labelled  selectively  as  shown 
in  Figures  4.10  and  4.11.  The  two  new  nodes  located  along  the  interface 
between  superelements  1  and  4,  namely  nodes  160  and  162,  remain  con¬ 
strained  as  do  nodes  167  and  146  located  between  superelements  2  and  3 
since  superelements  3  and  4  are  not  refined.  These  four  nodes  are  inter¬ 
face  nodes  and  are  marked  as  negative  in  the  restart  data  produced  under 
the  REFINE  header  in  the  final  output  file.  Chapter  5.0  contains  a 
description  of  the  reason  for  this  procedure. 

Nodes  145  and  161  of  superelement  1  and  nodes  166  and  147  of  super¬ 
element  2  are  identified  as  slave  nodes.  These  nodes  would  become  master 
nodes  if  boundary  conditions  were  applied  to  them.  The  remainder  of  the 
nodes  labelled  as  slave  nodes  are  internal  to  each  substructure  and  are 
constrained  since  their  dependent  basis  is  contained  in  another  element. 
The  new  master  nodes  which  are  identified  as  a  result  of  the  refinement 
are  the  new  nodes  located  along  the  interface  between  superelements  1  and 
2  and  also  the  four  existing  nodes  required  to  be  identified  as  master 
nodes  for  each  refined  superelement  because  these  nodes  were  identified 
in  the  dependent  basis  of  a  constrained  nodes. 


NGN , NBORD , NDDOF , IDCNK , XO , YO , ZO 


—  BORD  NO.,  NMNP,  (MMM(I) ,1=1 ,NMNP) 

NSS  NBORD 

where:  MMM  contains  master  node 
numbers  the  boundary  node 
1—  is  dependent  upon 

NSUB,NMN, NBORD 

(NMN1 (I) ,1=1 ,NMN) 

(NMN2(I) ,1=1 ,NMN) 

(SPACE(I), 1=1, NBORD) 

where:  NMN1  contains  substructure  node  numbers 

NMN2  contains  master  node  numbers  and  boundary 
nodes  are  input  as  negative  number 
SPACE  contains  boundary  node  numbers 
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TABLE  4.3 

Superelement  Data  for  Unrefined  Bulkhead  Model 
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Superelement  Data  for  Refined  Bulkhead  Model 
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62 

0. OCOE  +  OO  O.OOOE  +  OO  O.OOOE+OO 
O.OOOE+OO  0. OCCE+CO  O.OOOE+OO 
O.COOE+OO  0. OCOE+OO  O.OOCE+OO 
O.OOCE+OO  O.OOOE+OO  O.OOOE+CO 
-0.12SE+03  0. 125E+01  O.OOCE+OO 
0 . 1 25E+03-0 . 1 25E+0 1  O.OOOE+OO 
-C.250E+00  C.75CE-02  O.OOOE+OO 
-0.25CE+00  0.750E-02  O.OOOE+CC 
0 . 750E-02  0.500E+00  0. OCOE+OO 
0.75CE-02  0 . 500E+0C  O.OCCE+OC 
O.OOOE+OO  O.OOOE+OO  O.COOE+CO 
O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 

O.OOCE+OO  O.OOOE+OO'  0. OCOE+OC 
C.OOOE+OO  0. OCOE+OC  O.OOOE+CO 
0. OCOE+OO  O.CCGC+CC  O.OOOE+CC 
O.COOE+OO  0. OCOE+OC  0. CCOE+CC 
-0. 109E+C3-0. 125E+C1  C. OCOE+OC 
0. 109E+03  0. 125E+01  O.OOOE+CO 
-0 . 250C+C0-C . 862E-C2  O.OCCE+OO 
-0 . 250E+0C-0 . S52E-02  O.OOCE+OO 
-0 . 882E-02  0 . 5CCE  +  CC  O.OOCE+OO 
-0.8S2E-02  0 . SOOE+OC  O.OCCE+OO 
O.OCCE+OC  O.OCCE+OC  0. OCOE+CC 
O.OCCE+OC  0. COCE+OO  0. OOOE+OC 

O.OOOE+OO  O.OOCE+OO  C.OOOE+CC 
O.OCCE+OC  0. OCCE+CO  0. OOOE+OC 
O.OOOE+CO  O.OOCE+OO  0. OCOE+OO 
O.OOOE+OO  O.OOOE+OO  O.OOOE+OO 
0 . 1 09E+03  O.COOE+OO  0. OOOE+OC 
■0.109E+03  0. OOOE+OC  O.OCCE+OO 
■C.25CE+00  O.OOOE+CO  O.OOCE+OO 
0 . 25CE+0C  O.OCCE+OO  0. OCOE+OC 
O.OCOC+OO  0 . 50CE+CC  0. OOOE+OC 
O.OCCE+CC  C.5CCE+C0  C.CCCE+CC 
C.OOOE+CC  C.OCCC+CC  C.CCCE+CC 
G.CCOE+QO  0. OOOE+OC  O.OCCE+OC 

O.OGOE+OC  0. OCOE+OC  0. OCOE+OC 
O.OOOE+OO  O.OOOE+CO  O.COOE+CO 
C. OCOE+OC  C. OCCE+CO  C.OOOE+CC 
O.OOCE+OO  0. OCOE+OO  0. COCE+OO 
0. 10SE+03  C.125E+01  C.OOOE+OO 
C. 1C5E+C3-C. 125E+01  O.COOE+OO 
0 . 25CE+00  0.882E-02  O.OOOE+OO 
0 . 250E+00  C.382E-C2  O.CCCE+CO 
0 . 882E-02  0.5CCE+CC  O.CCCE+OC 
0 . 382E-02  0 . 500E+0C  O.OOOE+CO 
O.OOOE+OO  O.OOOE+OO  0. OCOE+OO 
O.OOOE+OO  O.COOE+OO  O.COOE+OO 


15  16 

31  32 

47  48 

1  C  7  11C 

15  16 

31  32 

47  48 

53  54 
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’29  2  5 

89  C.500E+00  O.OCOEi-Oa  O.OCOC+OO  0. OCOE+OO  0. OOCE  +  OO  0. COCE+CO 

94  0.500E+00  0. COCE+CO  C-OCCE+OO  O.OOOE+CO  O.OOOE+OO  O.OOOE+OO 

39  O.OOOE+OO  0.500E+00  O.OCCE+OC  O.OOOE+OO  O.COOE+OO  O.OOOE+OO 

94  O.OCOE+CO  0 . 50CE+C0  O.OOCE+OC  O.OOOE+OO  O.OOOE+CO  O.OOOE+OO 

39  O.OCOE+CO  O.OOCE+OC  0 . 500 E +00 -0 . 38 7 C+02  0-25CE+01  O.CCCE+OO 
94  O.OOOE+OO  0 . GOCE+CO  0.500E+00  O . 98 7E+02 -0 . 25CE+0 1  O.OOCE+OC 
83  O.COOE+OO  0 . OOOE+OC  0 . 1 9CE-02 -0 . 25CE+0C  C.19CE-01  O.OOOE+OO 
94  O.OOOE+CO  0 .OOOE+OO-O . 130E-02-0 . 25CE+00  0.190E-01  O.OOOE+CO 
39  0 . OCOE  +  OO  „C . OCOE +00-0 . 430E-C4  0.1900-01  C.SOCG  +  OO  O.OOOE  +  OO 
94  O.OCOE+CO  0 . COOE+CO  0-48GE-C4  0.190E-01  0.500E+00  O.COOE+OO 

39  C-CCCE+OC  O.OOCE+OC  O.OCOE+CO  C.OCOE+OO  C.OCOE+OC  O.OOOE+OO 

94  O.CCCE+OO  0 . COCE+CO  0 . OCCE+OO  O.OCOE+CO  C.OCOE+OO  O.OOOE+OO 

130  2  5 

12  0.50  OE  +  OG  O.OOCE  +  OC  O.OOCE  +  OC  C.OOOE  +  OO  C.COCE  +  OO  O.OOOE  +  OO 

94  0 . 50CE+00  0. COOE+CO  O.OOCE+OC  O.OOOE+OO  C.COCE+OO  O.OOOE+CO 

12  C.OOCC+OC  0 . 500E+00  O.OCOC+OO  C.OCOC+OO  O.OOOE+OO  0. OOOE+OC 

94  0. OCCE+OO  0 . 500E+00  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.COOE+OO 

12  O.OOOE+CO  0 . OCCE+OO  0.500E+00  0.101E+03  O.OCCE+CO  O.OOCE+OC 

94  O.OCOE  +  CO  O.OOOE  +  OO  0 . 5C0E  +  00-0 . 1 0 1 E  +  03  O.OCOE  +  OC  O.OCOE  +  CO 

12  0 . COOE+OC  0. OOOE+OO-O. 1 85E-02-0 . 25CE+00  O.OCOC+OO  C.OOCC+OC 
34  0 . COOE+OO  0 . OCCE+OO  0 . 1 85E-02-0 . 250E+C0  O.COOE+OO  O.OOOE+OO 
12  O.OOOE+OO  O.OOOE+OO  O.OCOE+CO  O.OOCE+OC  0.500E+00  O.OCOE+OC 

94  O.OOOE+OO  O.OOOE+OO  O.OCCE+OC  O.OOOE+OO  0.500E+00  O.OOCE+OC 

12  0. OCCE+OO  O.COOE+OO  C.CCOE+OO  O.OCOE+OC  O.COCE+OC  O.OOOE+OO 

34  O.COCE+OC  O.OCCE+OC  O.CCCE+OO  O.OOOE+OO  O.OOOE+CO  O.OOOE+OO 

145  2  8 

37  0 . 500C+00  O.CCCE+OO  O.OOOE+OO  O.OOOE+CO  O.OOOE+OO  O.OOOE+OO 

38  0 . 500E  +  00  C. OCCE+OO  0. COCE  +  CO  O.OOOE+OO  O.OOOE  +  OO  0. COOE  +  OC 

37  O.OOOE+OO  0.500E+00  O.OOOE+OO  O.OOOE+CO  O.OOOE+CO  O.OCOE+OC 

38  O.OOOE+OO  0.500E+00  O.OOOE+OO  O.GOOE+CO  O.OOOE+OO  O.OOOE+OO 

37  C.OOOE+OO  0 . OCOE+OO  0.500E+00  O.COOE+OO  0.70CE+02  0. OCCE+OO 

38  O.OOOE+OO  0 . OOCE+OG  0.50CE+00  C . OCCE+CO-O . 7 OC E+0 2  O.OOOE+OO 

37  O.OOOE+OO  0 . OOCE+OO  O.OOOE+OO  0.5CCE+0C  C.OOCE+CO  O.OOOE+OO 

38  O.OOOE+OO  O.OOOE+OO  C. OOCE+OO  0.500E+00  0. OCCE+OO  C.OOOE+OO 

37  O.OOOC+OO  0 . OOOE+OO-O . 268E-02  0 . OOOE+O 0-0 . 2 5CE+0 0  0. OCOE+OO 

38  O.OOOE+OO  O.COOE+OO  0.2S8E-02  0 . OOOE+CO - 0 . 2 50 E +00  O.OOOE+CO 

37  O.OOOC+OO  C.OOOE+OO  C.OCOC+OO  O.OOCE+OC  O.OCOE+OC  O.OOOC+OO 

38  0. OCOE+OO  O.COOE+OO  0. OCCE+OO  O.COOE+OO  O.COOE+OO  O.OOOE+CC 

ISO  2  6 

13  C.500E+00  0. OOCE+OO  0. OCOE+OO  O.OOOE+CO  0. OCCE+OO  O.OOOE+OO 

14  0 . 500E+C0  0. OOCE+OO  O.OOOE+OO  0. COCE+CO  O.OOOE+CO  O.OOOE+CO 

13  C.OOCE+CO  C.50CC+OC  O.OCOC+OO  O.OOCE+OC  0. OCCE+OO  O.OCOE+OC 

14  0. OOCE+OO  0.50CE+Q0  O.OOOE+OO  O.OOOE+OO  0. COOE+CO  C.OOOE+OO 

13  C.CCOE+OO  C-OCOE+OO  0.50CC+00  0.12SE+01  0.S74E+02  O.OCOE+CO 

14  O.OOOE+CO  C.OOCE+CO  0 . 5CCE+C0-0 . 1 25E+C 1 -0 . S74E+02  O.OCCE+OC 

13  0. OOCE+OO  0 . OCOE+OC-O . 5 1 SE-C4  C . 500C + 0 C -0 . 1 39 C-C 1  O.COCE+OC 

14  O.OCOE+CO  0 . OOCE+OO  0.51SE-C4  0 . 5C0E+C0-0 . 1 3 9E-C 1  O.OCOE+OC 

13  O.COOE+OO  O.OOCE+00-0.278E-02-0.139E-01-0.250E+00  O.OCOE+CO 

14  O.OOOE+OO  0. OCOE+OO  0 . 27 8E-02 -0 . 1 39E-C 1 -C . 2 5 CE+CC  O.OCOE+OC 

13  0. OCCE  +  OO  O.OOOE  +  CO  0. OCOE  +  OO  C.OOOE  +  OO  C.CGCE+.CC  O.OCOE  +  OC 

14  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOCE+CO  O.CCCE+OO 

161  2  6 

36  0  - 500E  +  00  O.OOOE  +  OO  O.OCOE  +  OC  C.OCOC  +  OO  0. OOCE  +  OO  0. OCCE  +  OO 

37  0.500E+C0  O.OOOE+OO  0. OOOE+OC  0. COOE+CO  O.OOOE+OO  0. OCOE+OO 

36  O.OOOE+OO  0.500E+00  0. OCOE+OO  O.OCOE+CO  O.OOCE+OC  O.OOOC+OO 

37  0. OCOE+OO  0.500E+00  O.OOOE+CO  C.OOOE+OO  O.OOOE+OO  O.OCOE+CO 

36  O.COOE+OO  O.OOCE+OC  0.500E+00  O.OOOE+OO  0.S68E+02  0. OCOE+OO 

37  0. OOCE+OO  0. OOCE+OO  0.50CE+00  0 . 000 E+ 0 C-0 . 6 6 8 E+02  O.OOOE+OO 

36  O.OOOE+OO  O.OOCE+OC  O.OOOE+OO  0.500E+00  O.OOOE+OO  O.OCOE+OC 

37  O.OOOE+OO  C. OOCE+OO  O.OCOE+OC  0.500E+00  0. OOCE+OO  O.OCCE+OC 

36  O.OOOE+OO  0 . OOOE+OO-O . 28 1 E-02  0 . OOOE+OO-O . 250E+00  O.OOOE+OO 

37  0. OOCE+OO  O.OOOE+OO  0.281E-02  0 . OOOE+OO-O . 250E+00  0. OOOE+OC 

36  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  0. OOCE+OO  0. OCOE+OO 

37  O.OOOE+OO  0. OOCE+OO  O.OCOE+OC  O.COOE+OO  O.OOOE+OO  0. OOOE+OC 

162  2  6 

12  0 . 500E+00  O.OCOE+OC  O.OOOC+OO  O.OOOE+OO  0. OCCE+OO  O.COOE+OO 

13  0 . 500E+00  O.COOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OOOE+CO 

12  O.OOOE+OO  0 . 500C+00  O.OOOE+OO  O.OOOE+OO  O.OOOE+OO  O.OCOE+OC 

13  0. OOOE+OC  0.500E+0C  0. OCOE+OO  O.OOOE+OO  O.OOCE+OC  O.OCOE+OC 

12  0. OOOE+OC  O.OOOE+OO  0.500E+00  O.OOOE+OO  0.691E+02  O.OOOE+OO 

13  C.OCOE+OO  O.OOOE+OO  0.5C0E+C0  0 . COOE+C 0 -0 . 69 1 E+02  C.OOOE+OO 

12  O.COOE+OO  O.OOOE+OO  O.CCCE+OO  0.5CCE+00  O.OOOE+OO  0. OCOE+OO 

13  O.OOOE+OO  O.OOOE+OO  0. OCOE+OO  C.500E+00  O.OOOE+OO  O.COOE+OO 

12  O.OOOE+OO  0. OOOE+OO-O. 271E-02  C . 0 OOE+O C- 0 . 250 E + 00  O.OOOE+OO 

13  O.OOOE+OO  O.OOOE+OO  0.271E-02  0 . COCE+OO-O . 25CE+00  O.OCOE+CO 

12  0. OOCE+OO  O.OOOE+OO  O.COOE+OO  O.COOE+OO  O.OCOC+OO  C.OOOE+OO 

13  O.OOOE  +  OO  O.OOOE+OO  0. OCOE  +  OO  O.OOOE  +  OO  O.OOCE  +  OC  O.OCOE  +  CO 


P151394.PDF  [Page:  42  of  80] 


TABLE  4.4  (Continued) 

r  •  V  ft- 


4.17 


2 

91 

10 

. 

1 

2 

3 

4 

5 

5 

7 

3 

9 

1  0 

1  7 

18 

1  9 

20 

21 

22 

23 

24 

25 

26 

33 

34 

35 

.  36 

37 

33 

39 

40 

4  1 

4  2 

49 

50 

51 

52 

53 

54 

S5 

56 

57 

69 

108 

107 

1  08 

103 

1  10 

1  1  1 

VI  2 

1  13 

114 

1  1  5 

67 

72 

89 

94 

144 

145 

160 

162 

1  63 

1  64 

65 

55 

67 

S3 

69 

70 

71 

72 

73 

74 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

96 

97 

98 

93 

100 

38 

101 

102 

103 

104 

1  1  1 

1  1  2 

51 

"1  1  3 

1  14 

1  1  5 

1  IS 

1  17 

1  18 

58 

1  22 

1  23 

124 

125 

126 

1  27 

1  28 

129 

1  30 

131 

278 

279 

280 

281 

287 

238 

285 

286 

282 

2.8  3 

0. OOOE+OO 

0  -  OOOE+OO 

0. OOOE+OO 

. 

124 

2 

6 

1 i 
27 
43 
7  0 
1  16 
1  65 
75 
30 
105 
59 
132 


1  2 
28 
44 
91 
1  17 

75 
31 
105 
SO 
1  33 


1  3 
23 
45 
92 
1  18 

77 

92 

107 

61 

134 


3 5  0.50CE+0C 
67  C.5C0E+CC 

36  0 . OOOE+OC 
67  Q.QCOE+OO 


67  0. OOOE+OO 
3.6  0  .  OOOE  +  OC 
S7 
125 


0. OOOE+OO  0. OOOE+OO 
0. OOOE  +  OO  0. OOOE  +  OO 
0 . 50GE  +  C0  0. OOOE  +  OO 
0 . 500E+CC  0. OOOE+OC 
36  0 . OOCG+OC  0. OOOE+OO  0.500E+00- 
67  0. OOOE+OO  0 . OCOE+OO  0.500E+0C 
36  0. OCOE  +  OO  0. OOOE+OO  C.150E-C2- 
67  0 . OCCE+OC  0 . OOCE+OO-O . 1 50E-02- 
36  O.OOOE+CO  0. OOOE+OC  0 . 1 50G-04- 
0  .  OOOE  +  OC-O 
0 .  OOOE  +  OO  0 

0. OOOE+OO  0. OOOE+OO  0. OOOE+OO 
2  5 

67  0.500E+00  0. OOOE+OO  0:0OCE+00 
72  0 . 500E+0C  0. OOOE+OO  O.OOCE+OO 
57  O.OOCE+OO  0 . 5C0E+00  0. OCOE+OO 
72  C. OOOE  +  OC  0 . 500E  +  0C  O'.OOOE+OO 
57  0 . OOOC+OC  0. OOOE+OC  0.500E+00- 
72  0. OOOE+OO  0. OOOE+OO  0.500E+00 
67  O.OCCE+OC  0. OOOE+OO  C.172E-C2- 
72  C. OCOE+OO  C.COOE+CO-O. 172E-Q2- 


.  1 50E-04- 
.  COOE  +  OC 


0. OOOE+OO  0 
0. OOOE+OO  0 
O.OOCE+OO  0 
0. OCOE+OO  0 
0 . 1 25E+03-0 
0 . 1 25E+03  0 
0 . 250E+00-0 
0 . 250E+CC-0 
0 . 750E-02  0 
0 . 750E-C2  0 
0. OOOE+OO  0 
0. OCOE+OO  0 


0 -OOOE+OO 
0. OOOE+OO 
0 .OOOE+OO 
0. OOOE+OO 
0 . 1 0SE+03 
0  -  1 09E  +  03-0 . 
C.25CE+0C  0. 
C.2S0E+0C  0. 


.  OCOE+OO 
.  OCOE  +  OO 
-OOOE+OO 
.  OOOE  +  OO 
.  1 2SE  +  01 
.  125E+01 
.  750E-02 
.  750E-02 
.  500E+CC 
.  500E  +  00 
•OOOE+OO 
.  OCOE  +  OO 

•OOOE+OO 
.  COOE  +  OO 
.  OOOE+OO 
OOOE+OO 
.  125E  +  01 
125E+01 
.  862E-02 
.  862E-G2 


Q.OCCE+OO 
0 . OOOE+OO 
0 -OOOE+OO 
0. OCOE+OO 
0. OCOE+OO 
0 . OCOE+OO 
0 . OCOE+OO 
O.OCCE+OC 
0. OOOE+OO 
0. OCOE+OO 
0  - OCCE  +  OO 
0 -OOOE+OO 

0. OCOE+OO 
0 . OOOE+OO 
0. OCOE+OO 
0 . OOOE+OO 
0 . OCOE+OO 
0 . OOOE+OO 
0 . OCOE+OO 
0 -CCOE+CO 


67 

72 

67 

72 

126 

49 

72 

49 

72 

49 

72 

49 

72 

49 

72 

49 

72 

130 

49 

89 

49 

89 

43 

83 

49 

83 

49 

89 

49 

89 


0. OOOE+OC  0. OOOE+OO 
0. OOOE+OO  0. OOOE+OO 
0. OOOE+OC  C.OCOE+CO 
C. OOOE+OO  0. OCOE+OO 
2  6 

C.500E+0C  0. OCOE+OO 
0 . 500E+00  O.OOCE+OO 
0. OCOE+OO  0.500E+00 
C . OOCE+OO  0 . 50CE+C0 
O.OOCE+OO  0. OCOE+OO 
0. OOOE+OO  O.OOCE+OO 
0. OOOE+OO  0. OOOE+OO 
O.OOCE+OO  O.OOOE+CO 
0. OOOE+OO  0. OCOE+OO 
0 . COOE+OG  0. OOOE+OO 
0. OOOE+OO  0. OOOE+OO 
0. OOOE+OO  0. COOE+OO 
2  S 


C.1S8E-04  C.362E-02 
0.198E-04  0.862  E- 0  2 
O.OCOC+OO  O.OOCE+OO 
0. OCOE+OO  0. OOOE+OO 

0. OOOE+OO  0. OCOE+OO 
0. OCOE+OO  O.OOOE+CO 
0. OCOE+OO  0. OOOE+OC 
0. OOOE+OO  0. OOOE+OO 
0.500E+00  0.1C9E+03 
0 . 50CE+00-0 . 1 03E+03 
0 . 1 72E-02-0 . 2S0E+00 
0. 172E-02-0.25CE+00 
C. OOOE+OO  C. OOOE+OO 
O.OOOE+CO  O.OOCE+OO 
0. OOOE+OO  0. OOOE+OO 
O.OOCE+OO  0. OOOE+OO 


0 . 500E+0C  O.OOCE+OO 
0 . 5C0E+CC  0. OOOE+OC 
0. OOOE+OO  0. OOOE+OO 
0 . OOCE+CC  0. OCOE+OO 

O.OOCE+OO  0. OOOE+OO 
0. OCOE+OO  C. OOOE+OO 
0. OOOE+OC  0. OOOE+OC 
O.OOOE+CO  0. OCOE+OO 
0. OCOE+OO  0. OOOE+OC 
0 . OCCE+OO  0. OOOE+OC 
0. OCOE+OO  0. OCOE+OO 
0. COOE+OO  O.OOCE+OO 
0 . 5C0E+0C  0. OOOE+OC 
0 . 500E+00  O.OOCE+OO 
0. OOOE+OC  0. OCOE+OO 
0. OOOE+OO  0. OCOE+OO 


0 . 500E+00  C. OOOE+OO  O.OOCE+OO  O.OOCE+OO  0. OCOE+OO  O.COOC+OO 
0 . 500E+00  0. OOOE+OO  0. OOOE+OO  0. OOOE+OO  O.CCOE+OC  0. OCCE+OO 
0. OOOE+OO  0 . 5C0E+00  0. OOOE+OO  0. OOOE+OO  0. OOOE+OO  0. OOOE+OO 
O.OOCE+OO  0.500E+00  0. OOOE+OO  0. OOOE+OO  0. OOOE+OC  0. OOOE+OO 
C.OOCE'+OO  0. OOOE+OO  0 . 5  0  0  E+ C  0 -0 . 1  0  6  E  + 03 -0 . 1  2  5E+ C  1  0. OOOE  +  OO 
0. OOOE+OO  0. COOE+OO  0.500E+00  0.10SE+03  0.125E+01  0. OOOE+OO 
O.OOCE  +  OO  0 -  OOOE  +  OO  0 . 1 7 SE- 02 -0 . 25 0 E  + 00 -0 . 3 8 2E- 02  0. OCOE  +  OO 
0. OOOE+OO  0 . OOCE+OO-O . 1 76E-02-0 . 25CE+00-0 . 882E-02  0. OOOE+OO 
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FIGURE  4.2:  Flowchart  for  Program  REFINE 
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FIGURE  4.2:  Flowchart  for  Program  REFINE  (Continued) 
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:  Flowchart  for  Program  REFINE  (Continued) 
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FIGURE  4.4:  Substructured  Bulkhead  Model  Before  Refinement 
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FIGURE  4.5:  Substructured  Bulkhead  Model  After  Refinement 
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FIGURE  4.6:  Master  Nodes  for  First  Superelement 
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FIGURE  4.7:  Master  Nodes  for  Second  Superelement 
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FIGURE  4.8:  Nodes  for  First  Superelement  After  Refinement 
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FIGURE  4.9:  Nodes  for  Second  Superelement  After  Refinement 
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FIGURE  4.10:  Slave  Nodes  for  First  Superelement  Defined  Interactively 
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FIGURE  4.11:  Slave  Nodes  for  Second  Superelement  Defined  Interactively 
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5.1 


CHAPTER  5 

RESTART  CAPABILITY 


5.1  Introduction 


The  REFINE  program  may  eventually  be  automated  to  the  point  that  it 
could  operate  as  a  component  of  an  adaptive  finite  element  model 

optimization  capability  which  once  initiated  would,  without  any  user 
intervention,  alternately  refine  the  mesh  and  resolve  until  desired 
levels  of  accuracy  are  attained  with  near  optimal  distribution  of 
freedoms  in  the  finite  element  model.  But  considerable  further 
development  of  the  algorithms  within  both  the  discretization  error 

estimation  and  mesh  refinement  capabilities,  is  still  required.  The 

capability  to  refine  a  previously  refined  structure  provides  a  so  called 
"restart"  feature  and  moves  a  step  closer  to  the  completely  automated 
procedure . 

Prior  to  the  development  of  this  capability,  the  user  had  to 

reconsider  the  original  model  and  refine  the  model  to  increasing  levels 
as  desired.  Now  the  user  has  only  to  "restart"  and  refine  the  previously 
refined  model  in  areas  of  interest.  Considerable  savings  in  time  are 
possible  since  the  user  no  longer  has  to  reconstruct  the  refinement 
process  to  attain  the  starting  point  from  which  to  perform  the  desired 
refinements.  The  implementation  of  this  capability  is  discussed  in  this 
chapter . 

5.2  Implementation 

The  refinement  of  an  element  involves: 

1 )  generation  of  new  nodes ; 

2)  replacement  of  the  connectivity  of  a  parent  element  with  the  connec¬ 
tivities  of  elements  produced  by  the  refinement;  and 

3)  generation  and  manipulation  of  any  resulting  constraint  equations. 
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When  restarting  the  refinement  of  a  finite  element  model,  no  particular 
problems  are  expected  with  the  first  two  steps  but  difficulties  may  arise 
with  the  manipulation  of  constraint  equations.  The  algorithm  which 
determines  whether  constraint  equations  remaining  on  each  new  node 
requires  that  constraint  equations  be  available  in  the  uncondensed  form. 
This  implies  that  the  constraint  equations  for  the  original  model  must  be 
available  in  uncondensed  form.  The  algorithm  for  evaluating  whether 
constraint  equations  remain  on  a  node  is  explained  in  detail  in  Chapter 

4. 


The  uncondensed  multipoint  constraint  data  for  a  refined  model  can¬ 
not  be  readily  extracted  from  the  condensed  constraint  equations  data 
written  to  the  stiffness  modification  data  section  of  the  input  file. 
The  key  to  restarting  the  refinement  of  a  finite  element  model  is  just 
simply  to  store  the  uncondensed  constraint  equations  for  any  refined 
model  in  case  the  refinement  process  is  to  be  restarted.  The  restart 
data  associated  with  the  refined  model  data  is  therefore  stored  on  the 
USE  file  under  the  header  REFINE.  On  restarts,  the  REFINE  program 
searches  for  this  header  and  transfers  the  data  to  a  scratch  file  prior 
to  any  constraint  data  processing. 

Additional  data  is  required  when  restarting  the  refinement  of  a  sub¬ 
structured  model  if  the  nodes  which  are  boundary  nodes  are  to  be  identi¬ 
fied.  Recall  that  a  boundary  node  is  one  which  is  located  on  the 
perimeter  of  a  superelement.  Boundary  nodes  are  flagged  under  the  REFINE 
header  by  a  change  of  sign  on  the  dependent  node.  A  new  subroutine 
RESTART  was  developed  to  extract  the  information  about  the  boundary  nodes 
from  the  ''restart"  data  and  to  fill  the  master  node  and  boundary  node 
arrays.  Master  node  and  boundary  node  data  are  used  in  subroutine  MPCSE 
to  determine  whether  a  constraint  equation  on  a  node  is  to  be  retained 
and  the  node  becomes  a  slave  node  or  alternately  the  constraint  equations 
on  a  node  are  eliminated  and  the  node  becomes  a  master  node. 
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The  format  of  the  data  under  the  REFINE  header  is  the  same  for  each 
superelement  of  a  substructured  model  as  for  an  unsubstructured  model. 
The  variable  NDDOF  equal  to  the  number  of  degrees  of  freedom  is  entered 
as  zero  for  those  superelements  which  do  not  contain  any  slave  nodes. 

5.3  Sample  Problem 

The  sequence  of  steps  involved  in  the  recursive  refinement  of  a 
simple  model  composed  of  only  two  four-noded  shell  elements  is  now 
described.  This  example  illustrates  the  difference  between  "uncondensed " 
and  "condensed"  equations  and  the  requirement  for  the  former  in  the 
restarted  refinement  of  the  model.  The  orginal  model  composed  of  two 
four-noded  shell  elements  is  shown  in  Figure  5.1(a). 

The  first  step  in  the  refinement  of  the  model  is  to  refine  one  ele¬ 
ment,  as  shown  in  Figure  5.1(b),  and  generate  a  constraint  equation  for 

node  c  in  subroutine  REFINE  with  nodes  a  and  b  specified  as  dependent 
nodes . 


One  of  the  new  elements  is  then  refined  further  as  shown  in  Figure 
5.1(c).  This  results  in  the  creation  of  constraint  equations  for  nodes 
e,f,  and  g.  Node  e  is  dependent  upon  nodes  a  and  c,  node  f  is  dependent 
upon  nodes  c  and  d,  and  node  g  is  dependent  upon  nodes  d  and  h.  Because 
node  c  is  already  a  dependent  node,  the  constraint  equations  for  nodes  e 
and  f  are  said  to  be  "nested". 

After  all  specified  elements  have  been  refined,  the  final  step 
dealing  with  constraint  equations  involves  determining  which  constraints 
are  to  be  retained  and  writing  them  to  the  final  output  file  in  a  format 
acceptable  to  VAST.  In  the  example  being  considered,  the  constraints  on 
nodes  c,f,  and  g  will  remain  since  independent  basis  of  each  is  contained 
in  connectivity  of  another  element.  The  constraint  on  node  e  also 
remains  since  it  is  dependent  upon  a  dependent  node  and  is  not  an  inter— 
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face  node.  The  reader  will  recall  that  by  the  convention  established  in 
the  introduction  to  the  present  chapter,  an  interface  node  is  defined  as 
a  node  located  on  the  perimeter  of  a  parent  element  that  is  common  to  two 
or  more  elements.  The  constraint  equation  for  nodes  e  and  f  are  conden¬ 
sed  so  that  no  dependent  nodes  are  involved.  The  final  result  is  node  e 
dependent  upon  a  and  b  and  node  f  dependent  upon  nodes  a,  b  and  d. 

Figure  5.1(d)  shows  the  result  of  a  restart  and  in  which  the  second 
of  the  original  elements  is  refined.  Before  the  restart  capability  was 
implemented  in  the  REFINE  program,  the  constraint  on  e  and  f  would  both 
be  incorrectly  removed  when  the  second  element  was  refined.  The  reason 
for  this  is  that  the  constraint  equations  on  node  c  would  be  removed  in 
subroutine  REFINE  since  the  independent  basis  for  this  equation  would  be 
determined  to  be  only  in  the  element  being  refined.  The  constraints  in 
the  remaining  nodes  would  again  be  investigated  by  subroutine  MPC.  The 
constraint  equations  on  node  g  would  remain  since  its  basis,  as  in  the 
first  model  refinement,  is  contained  the  connectivity  of  an  element.  The 
contraint  equations  on  nodes  e  and  f  are  however  removed  since  their 
independent  condensed  basis  are  not  in  the  connectivity  of  any  element  in 
the  model  and  there  are  no  dependent  nodes  in  the  equations. 


When  REFINE  is  restarted  using  the  uncondensed  equations,  the 
constraint  equations  are  manipulated  correctly.  The  constraints  on  nodes 
e  and  f  are  retained  as  they  should  since  independent  nodal  basis  is 
determined  to  be  in  an  element  connectivity. 
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(a)  Original  mesh  with  two  four-noded  shell  elements. 

a 


b 


(b)  Mesh  after  refining  one  of  the  original  elements. 
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(c)  Mesh  after  refining  one  of  the  elements  generated 
by  refinement. 


(d)  Mesh  after  refining  second  of  the  original  elements. 


FIGURE  5.1:  Recursive  Mesh  Refinements  and  Implication 
for  Dependent  Node  Determinations 
( O-  Dependent  Node ) 
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CHAPTER  6 

OTHER  EXTENSIONS  OF  SUPERELEMENT  REFINEMENT  CAPABILITIES 

6.1  Introduction 


Earlier  versions  of  the  program  REFINE  were  not  fully  operational  on 
superelements.  A  number  of  restrictions  had  to  be  adopted  to  facilitate 
producing  a  preliminary  capability  to  refine  superelements.  They  were 
identified  in  Chapter  4.  Items  d  and  e  were  addressed  in  Chapters  4  and 
5,  respectively.  Items  a  to  c  will  be  discussed  in  this  chapter. 

a.  Only  level  one  superelements  can  be  refined. 

b.  A  substructure  may  only  define  one  superelement. 

c.  Substructures  must  be  refined  in  ascending  order. 

6.2  Extension  to  Multilevel  Superelements 


The  restriction  that  only  level  one  superelements  can  be  refined 
still  applies.  The  extension  of  the  program  to  multilevel  superelements 
would  generalize  the  range  of  program  application  but  this  development 
appears  to  be  premature  at  the  moment.  It  is  preferred  that  the  super¬ 
element  capability  be  thoroughly  tested  and  evaluated  before  the  multi¬ 
level  superelement  capability  be  introduced. 

6.3  Removal  of  Limitation  that  One  Substructure 
Defines  One  Superelement _ 

The  removal  of  the  restriction  that  one  substructure  may  define  only 
one  superelement  represents  the  first  step  in  generalizing  the  refinement 
algorithm  for  superelements.  Although  necessary  program  modifications 
are  sketched  in  the  following  section,  they  have  not  as  yet  been  imple— 
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mented.  The  time  and  fiscal  constraints  of  the  current  contract  did  not 
permit  this  work  to  be  completed. 

Of  several  possible  approaches  to  perform  this  task  that  have  been 
investigated,  only  the  two  approaches  which  appear  most  feasible  are 
discussed  in  this  section. 

One  requirement  arising  from  the  removal  of  the  limitation  that  one 
substructure  define  one  superelement  is  the  need  for  an  accounting 
procedure  to  determine  the  number  of  new  substructures  created  to  relate 
superelements  with  them.  The  only  information  required  to  determine  the 
number  of  substructures  in  a  refined  model  is  a  list  of  superelement 
numbers  to  be  refined.  Prompting  the  user  at  the  start  of  the  program's 
execution  to  enter  the  superelements  to  be  refined  provides  data  which  if 
placed  in  an  array  can  be  used  to  direct  the  refinement  of  the  model. 

This  data  is  required  by  subroutine  RELEMB  which  reads  the  super¬ 
element  data  and  now  must  relate  the  superelement  to  the  correct  sub¬ 
structure  numbers.  The  VAST  data  input  requirement  that  superelement 
must  be  ordered  for  increasing  substructure  number  limits  the  options 
available  for  programming  the  superelement  refinement.  Firstly,  the 
original  superelement  number  scheme  can  remain  unchanged  and  the  sub¬ 
structure  numbering  changed  to  effectively  create  one  substructure  for 
each  superelement.  This  has  the  advantage  of  allowing  the  load  refine¬ 
ment  process  to  remain  unchanged  and  the  refinement  to  proceed  along  the 
input  file.  Secondly,  the  superelement  numbering  may  change  to  exploit 
any  repeated  substructures  which  remain  after  refinement. 

The  two  options  for  removal  of  the  limitation  that  one  substructure 
defines  one  superelement  are  best  presented  by  an  example.  The  example 
to  be  considered  involves  a  model  comprising  of  6  superelements  generated 
from  2  substructures  for  which  there  is  a  requirement  to  refine  super¬ 
elements  2  and  4  (see  Table  6.1).  The  new  substructure  numbering  scheme 
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resulting  from  changes  in  the  substructure  numbers  is  shown  in  Table 
6.2.  This  procedure  would  result  in  6  substructures  when  4  would  be 
sufficient  since  substructures  1,  3,  and  5  are  the  same.  The  renumbering 
of  the  superelements  only  results  in  the  minimum  of  4  substructures  as 
shown  in  Table  6.3.  However,  this  procedure  would  require  reordering  the 
load  data  as  well  as  the  superelement  data.  The  geometry  data  for  all 
the  substructures  to  be  refined  would  have  to  be  written  to  a  file  so 
that  the  refinement  process  could  proceed  sequentially  along  the  file. 

6.4  Refinement  of  Superelements  in  Arbitrary  Orders 

The  restriction  that  the  substructures  must  be  refined  in  ascending 
order  also  still  applies.  The  present  algorithm  proceeds  along  the  input 
file  with  no  rewinds  which  is  necessary  to  keep  the  I/O  at  a  minimum. 
The  execution  time  of  program  REFINE  appears  to  be  greatly  affected  by 
I/O. 
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TABLE  6.1 

Sample  Model  with  Six  Superelements 
Generated  from  Two  Substructures 


Super element 
Number 

Substructure 

Number 

1 

1 

2* 

1 

3 

1 

4* 

1 

5 

1 

6 

2 

*  Superelement  to  be  refined. 


TABLE  6.2 

Superelement  Refinement  with  all  Superelements 
Associated  with  Independent  Substructures - 


Superelement 

Substructure 

Number 

Number 

Old 

New 

Old 

New 

1 

1 

1 

1 

2* 

2 

1 

2 

3 

3 

1 

3 

4* 

4 

1 

4 

5 

5 

1 

5 

6 

6 

2 

6 

_ 

*  Superelement  to  be  refined 


TABLE  6.3 

Superelement  Refinement  with  Substructure  Reordering 
to  Minimize  Substructures  xn  Model - 


Superelement 

Number 

Old  New 

Substructure 

Number 

Old  New 

1 

1 

1 

1 

3 

2 

1 

1 

5 

3 

1 

1 

6 

4 

1 

2 

2* 

5 

1 

3 

4* 

6 

2 

4 

*  Superelement  to  be  refined 
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CHAPTER  7 

ADDITIONAL  DEVELOPMENTS  IN  PROGRAM  REFINE 


7.1  Introduction 


A  number  of  additional  changes  to  program  REFINE,  not  explicitly- 
called  for  in  the  contract  statement  of  work  but  determined  to  be  desir¬ 
able  during  development  work  and  in  application  of  the  program,  were  also 
implemented.  These  changes  are  described  briefly  in  the  following 
sections. 

7.2  Option  for  Data  Input  from  Separate  Files 

The  REFINE  program  was  modified  to  optionally  accept  the  following 
data  on  separate  files: 

1)  Geometry  Data  from  file  PREFX.GOM 

2)  Superelement  Data  from  file  PREFX.SED 

3)  Stiffness  Modification  Data  from  file  PREFX.SMD 

4)  Load  Data  from  file  PREFX.LOD 

The  REFINE  program  still  places  all  the  output  data  in  the  PREFN.USE  file 
and  not  in  separate  files. 

7-3  Arrays  Dimensioning  in  the  Main  Module 

The  changes  required  to  achieve  the  option  for  data  input  from 
separate  files  in  the  REFINE  program  resulted  in  some  more  general  organ¬ 
izational  changes  in  the  program.  For  instance,  array  data  transfers 
within  the  program  were  modified  so  that  all  arrays  are  now  dimensioned 
in  the  main  program  and  passed  through  the  argument  lists  of  the  various 
subroutines.  This  change  facilitates  any  future  increases  in  array 
dimensions  to  accommodate  large  scale  problems.  The  restrictions  adopted 
for  array  sizes  are  as  follows: 


P151394.PDF  [Page:  69  of  80] 


7.2 


NGNMAX 

=  999 

(number  of  geometric  nodes) 

NMNMAX 

=  999 

(number  of  master  nodes) 

NTIMAX 

=  999 

(number  of  time  steps) 

NERMAX 

=  250* 

(number  of  elements  to  be  refined) 

NEGMAX 

=  250* 

(number  of  element  groups) 

NDDMAX 

=  999 

(number  of  dependent  nodes) 

NIDMAX 

=  8* 

(number  of  independent  nodes) 

NDFMAX 

=  6* 

(number  of  dependent  degrees  of  freedom) 

NSKCMAX 

=  100 

(number  of  skewed  systems) 

All  of  the  above  array  sizes  except  those  marked  with  an  asterisk  corres 
pond  with  VAST  program  limitations. 

7.4  Superelement  Load  Data  Format  Change 

The  REFINE  program  was  modified  to  accommodate  the  new  format  in 
VAST  for  load  data  on  superelements  which  involves  the  card  identifying 
the  number  of  the  superelement  to  be  loaded  being  preceeded  by  an  aster¬ 
isk  (*).  The  implied  changes  for  the  REFINE  program  were  concentrated  in 
the  subroutine  RUSE,  which  reads  from  the  specified  input  file  and 
optionally  writes  it  to  the  specified  output  file  according  to  the  value 
of  the  ISW  parameter.  This  subroutine  previously  terminated  the  data 
transfer  process  and  returned  control  to  the  main  program  when  either  a 
"I" ,  "R”  or  was  found  in  the  first  column.  The  asterisk  previously 

only  identified  the  beginning  of  the  geometry  for  a  new  substructure. 
The  ISW  parameter  was  modified  so  that  "0"  indicates  that  the  data  is  to 
be  read  only,  "1"  indicates  that  the  data  is  to  be  read  and  written,  and 
"3"  indicates  that  the  asterisk  is  to  be  ignored  as  a  return  parameter. 


7.5  Identification  of  Nodes  Requiring  Multipoint  Constraint  Equations 

The  procedure  previously  employed  by  REFINE  to  determine  whether  a 
newly  generated  nodal  point  requires  a  constraint  or  alternately,  whether 
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the  constraint  on  a  node  can  be  removed  was  described  in  detail  in 
Appendix  E  of  Ref erence  1 .  The  only  check  to  be  performed  by  subroutine 
MPC  had  been  to  determine  if  any  independent  nodes  of  a  constraint 
equation  were  themselves  dependent  nodes.  If  they  were,  then  the 
constraint  equation  was  assumed  to  remain.  This  could  result  in  the 
retention  of  unnecessary  constraints  as  will  be  demonstrated  in  an 
example  problem  below. 

An  unrefined  model  for  a  square  plate  composed  of  18  VAST  triangular 
plate  bending  elements  (IEC=4)  is  considered.  The  central  element  is 
first  refined  to  order  2  and  then  the  resulting  four  elements  are  refined 
again  to  order  2.  The  final  mesh  is  shown  in  Figure  7-1.  The  refinement 
of  the  first  element  results  in  new  nodes  17,  18  and  19  which  are  all 
border  nodes  and  have  constraints  generated  for  them.  The  refinement  of 
the  first  of  the  new  elements  (connectivity  of  6,  17,  18)  results  in  the 
new  border  nodes  20,  21  and  22  with  constraints  generated  for  them. 
Similarly,  the  refinement  of  the  remaining  new  elements  produces  the  new 
border  nodes  23  to  28.  It  is  clear  that  although  the  constraint  equa¬ 
tions  on  nodes  22,  23  and  24  would  be  retained  as  a  result  of  their 
independent  nodes  being  themselves  dependent,  the  constraint  equations 
are  in  fact  not  necessary. 

As  explained  in  detail  in  Chapter  4,  this  algorithmic  deficiency  has 
now  been  corrected.  The  REFINE  program  was  modified  to  identify  and 
remove  the  constraint  on  these  nodes  by  using  the  concept  of  an  "inter¬ 
face"  node.  An  interface  node  has  been  defined  to  be  a  node  whose 
constraint  equation  is  common  to  two  or  more  parent  elements.  The 
additional  qualification  on  the  check  for  constraint  equation  retention 
is  that  if  dependent  nodes  are  found  in  the  constraint  equation  and  the 
node  is  not  an  interface  node,  the  constraint  should  remain. 

Returning  to  our  example  problem,  it  can  be  readily  demonstrated 
that  the  checks  for  constraint  equation  retention  is  now  reliable.  As 
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before,  the  refinement  of  the  centre  element  produces  new  nodes  17,  18 
and  19  which  as  border  nodes  have  constraint  equation  generated  for 
them.  Again  the  refinement  of  the  first  of  the  new  elements  produces  new 
border  nodes  20,  21  and  22  with  appropriate  constraint  equations  gener¬ 
ated.  The  refinement  of  the  second  new  element  results  in  new  border 
nodes  23  and  24  having  constraints  generated.  Node  22,  however,  is  spec¬ 
ified  as  an  "interface"  node  since  its  independent  basis  is  contained  in 
more  than  one  element.  The  refinement  of  the  third  and  fourth  new 
elements  results  in  new  border  nodes  25,  26,  27  and  28  having  constraints 
generated  and  23  and  24  being  flagged  as  "interface"  nodes.  (It  is 
important  to  note  that  the  logic  for  recursive  refinement  of  finite 
element  models  is  reliable  only  if  no  restarts  are  involved.  Sufficient 
information  is  not  stored  for  these , operations  to  be  carried  out  reliably 
if  model  refinement  is  restarted.) 

The  element  refinement  specifications  are  complete  and  now  sub¬ 
routine  MPC  is  used  to  determine  which  constraints  are  required. 
Constrained  nodes  17,  18  and  19  remain  since  the  comparison  against  the 
connectivity  list  of  the  independent  basis  determine  it  to  be  found. 
Constrained  nodes  20,  21,  25,  26,  27  and  28  remain  since  the  independent 
basis  of  these  nodes  contain  dependent  nodes  and  these  nodes  are  not 
interface  nodes.  The  constraints  on  nodes  22,  23  and  24  were  removed 
since  the  nodes  are  identified  as  interface  nodes  and  have  dependent 
nodes  in  their  constraint  equations. 

7.6  Graphical  Boundary  Condition  Selection 


When  program  REFINE  is  used  to  refine  a  finite  element  model,  new 
nodes  are  inserted  between  the  old  ones  and  any  regularity  of  node 
numbering  is  destroyed.  Without  extensive  experience  in  the  operation  of 
the  program,  the  user  will  not  be  able  to  readily  establish  node  numbers 
on  selected  portions  of  the  structure  as  required  for  instance  when 
defining  boundary  conditions.  Graphical  aids  of  some  sort  are  therefore 
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desirable,  the  minimum  being  an  option  to  plot  the  refined  model  or  a 
portion  of  it  and  to  label  the  node  numbers.  A  capability  which  is  some¬ 
what  more  sophisticated  and  automated  has  been  developed  within  REFINE  in 
a  previous  contract  [1].  It  utilizes  the  graphics  cursor  capabilities  of 
the  VASTG  graphics  program  PL0TV1 . 

The  graphical  node  point  selection  capability  is  optionally  acti¬ 
vated  when  defining  boundary  conditions  after  specifications  for  the 
refinement  of  a  substructured  or  unsubstructured  finite  element  model 
have  been  completed.  Within  the  current  contract,  the  subroutine  BCOND 
in  REFINE  program  has  been  extended  to  retain  boundary  conditions 
assigned  to  the  original  model  and  in  the  case  of  a  substructured  model 
to  ensure  that  a  node  having  boundary  conditions  applied  becomes  a  master 
node  (if  it  is  not  already  a  master  node). 

The  operation  of  the  program  to  define  problem  boundary  conditions 
is  described  in  detail  in  Appendix  D  of  reference  [1].  Basically,  the 
REFINE  program  transfers  control  to  the  VASTG  graphics  program  PLOT VI  to 
activate  graphical  node  selection.  This  selection  is  performed  by  plac¬ 
ing  a  "window"  around  the  nodes  of  interest.  However,  if  the  boundary 
conditions  are  to  be  specified  along  a  curved  edge  then  this  becomes  a 
tedious  procedure.  A  user  may  therefore  define  the  boundary  conditions 
using  the  VASTBC  [2]  program  after  the  model  has  been  refined.  VASTBC 
provides  the  user  with  the  option  to  select  nodal  points  individually 
using  the  graphics  cursor  or  multiply  using  the  graphics  cursor  to  define 
a  window  in  much  the  same  way  that  the  capability  in  REFINE  based  on  the 
program  PLOT VI  does. 

7*7  Treatment  of  Skewed  Coordinate  Systems 

When  boundary  conditions  are  associated  with  spatial  directions  not 
aligned  with  global  coordinate  directions,  a  skewed  coordinate  system  can 
be  utilized  to  align  model  degrees  of  freedom  with  directions  appropriate 
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for  the  definition  of  constraints.  Such  a  capability  has  been  implement¬ 
ed  within  VAST  [3]  and  thus  it  is  desirable  that  the  effect  of  defining 
skewed  coordinate  systems  be  completely  taken  into  effect  when  refining  a 
model.  The  REFINE  program  previously  had  the  capability  to  only  transfer 
the  skewed  coordinate  system  data  from  the  original  model  to  the  refined 
model.  In  this  contract,  REFINE  was  further  enhanced  to  account  for  the 
effect  of  the  skewed  coordinate  system  in  the  constraint  equation  defini¬ 
tions. 

7.8  Improved  Element  Selection  Capability 

The  REFINE  program  has  been  modified  to  provide  the  user  the  option 
to  select  elements  graphically  or  through  interactive  prompting  at  each 
stage  of  the  element  refinement  specifications.  The  program  previously 
allowed  the  user  only  once  to  define  the  option  for  the  selection  of  ele¬ 
ments  for  refinement  either  by  graphical  means  or  through  interactive 
prompting.  After  the  option  was  defined,  all  subsequent  element  selec¬ 
tions  would  have  to  utilize  the  same  selection  option. 

The  capability  to  select  elements  graphically  is  described  in 
Appendix  D  of  Reference  1 .  A  file  PREFX.ELM  with  the  elements  sorted  in 
ascending  order  with  respect  to  first  the  superelement  number,  next  the 
element  group  number,  and  finally  the  individual  element  number,  is 
produced.  The  organization  of  data  in  the  PREFX.ELM  file  in  ascending 
order  proved  extremely  beneficial.  In  REFINE,  it  was  possible  to  improve 
the  program  efficiency  (to  decrease  CPU  requirements)  and  to  reduce  the 
amount  of  I/O.  This  improvement  is  only  realized  in  the  cases  where  the 
window  contains  several  element  groups.  The  I/O  reductions  (and  associ¬ 
ated  time  savings)  were  accomplished  by  skipping  two  sections  of  the 
REFINE  main  program;  firstly,  the  one  that  copies  remaining  element  group 
data  to  the  output  file  and  secondly,  the  section  that  reverses  the  input 
and  output  file  unit  numbers  before  proceeding  to  refine  another  element 
The  refinement  proceeds  sequentially  through  elements  in  the 


group . 
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input  file  as  directed  by  the  PREFX.ELM  file  until  all  elements  are 
refined.  The  two  steps  which  were  skipped  were  generally  necessary  when 
selecting  elements  via  prompting  because  each  time  the  user  was  prompted 
for  an  element  group  to  be  refined,  the  reply  could  be  an  arbitrary  group 
number  and  could  therefore  involve  data  at  an  arbitrary  point  on  the  USE 
file. 


The  procedure  now  employed  to  specify  the  elements  through  inter¬ 
active  prompting  results  in  the  generation  of  a  PREFX.ELM  file  which  is 
the  same  format  as  the  file  created  when  the  elements  are  selected 
graphically.  The  refinement  proceeds  sequentially  through  the  elements 
identified  on  the  PREFX.ELM  file  and  therefore  the  unnecessary  data 
handling  is  eliminated. 


PRESSURIZED  PLATE  PROBLEM 
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FIGURE  7.1:  Recursively  Refined  Square  Plate  Model  (^=  Dependent  Modes) 
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